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Abstract 

This article reviews and updates the Standard Model prediction of the muon g—2. QED, 
electroweak and hadronic contributions are presented, and open questions discussed. 
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1 Introduction 



The evaluation of the Standard Model (SM) prediction for the anomalous magnetic 
moment of the muon = (gv — 2)/2 has occupied many physicists for over fifty 
years. Schwinger's 1948 calculation [1] of its leading contribution, equal to the one 
of the electron, was one of the very first results of QED, and its agreement with the 
experimental value of the anomalous magnetic moment of the electron, a e , provided 
one of the early confirmations of this theory. 

While a e is rather insensitive to strong and weak interactions, hence providing a 
stringent test of QED and leading to the most precise determination to date of the fine- 
structure constant a, a M allows to test the entire SM, as each of its sectors contribute 
in a significant way to the total prediction. Compared with a e , is also much better 
suited to unveil or constrain "new physics" effects. For a lepton I, their contribution to 
a>i is generally proportional to mf/A 2 , where mi is the mass of the lepton and A is the 
scale of "new physics", thus leading to an {m^/m e ) 2 ~ 4 x 10 4 relative enhancement 
of the sensitivity of the muon versus the electron anomalous magnetic moment. The 
anomalous magnetic moment of the r would thus offer the best opportunity to detect 
"new physics", but the very short lifetime of this lepton makes such a measurement 
very difficult at the moment. 

In a sequence of increasingly more precise measurements [2-5], the E821 Collab- 
oration at the Brookhaven Alternating Gradient Synchrotron has reached a fabulous 
relative precision of 0.5 parts per million (ppm) in the determination of a^, providing 
a very stringent test of the SM. Even a tiny statistically significant discrepancy from 
the SM prediction could be the harbinger for "new physics" [6]. 

Several excellent reviews exist on the topic presented here. Among them, I refer 
the interested reader to refs. [7-15]. In this article I will provide an update and a 
review of the theoretical prediction for a M in the SM, analyzing in detail the three 
contributions into which a^ M is usually split: QED, electroweak (EW) and hadronic. 
They are respectively discussed in sees. 2, 3 and 4. A numerical re-evaluation of the two- 
and three-loop QED contributions employing recently updated values for the lepton 
masses is presented in sees. 2.2 and 2.3. Comparisons between a'^ M results and the 
current experimental determination a^ xp are given in sec. 5. Conclusions are drawn in 
sec. 6. 

2 The QED Contribution to 

The QED contribution to the anomalous magnetic moment of the muon is defined as the 
contribution arising from the subset of SM diagrams containing only leptons (e, fj,, r) 
and photons. As a dimensionless quantity, it can be cast in the following general 
form [7, 16] 

a QED = Ai + A 2 {m ll /m e ) + A 2 (m fi /m T ) + A 3 (m M /m e , m M /m T ), (1) 

where m e , m^ and m T are the masses of the electron, muon and tau respectively. The 
term Ai, arising from diagrams containing only photons and muons, is mass indepen- 
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dent (and is therefore the same for the QED contribution to the anomalous magnetic 
moment of all three charged leptons). In contrast, the terms A 2 and A3 are functions 
of the indicated mass ratios, and are generated by graphs containing also electrons and 
taus. The renormalizability of QED guarantees that the functions A4 (i = 1,2,3) can 
be expanded as power series in a /it and computed order- by-order 

2.1 One-loop Contribution 

Only one diagram, shown in fig. 1, is involved in the evaluation of the lowest-order 
contribution (second-order in the electric charge); it provides the famous result by 
Schwinger [1], Af ] = 1/2 {A { 2 2) = A {2) = 0). 




Figure 1: Lowest-order QED contribution to a M . 
2.2 Two-loop Contribution 

At fourth order, seven diagrams contribute to a\ , one to A^\m^/m e ) and one to 
^2 (?rtu/m r ). They are depicted in fig. 2. The coefficient A± has been known for 
almost fifty years [17,18]: 

(4) = l^I + ^! + 3 c(3) _vr^ ln2 = _ Q 32g 478 9g5 579 / 3 n 
1 144 12 4 SU 2 ' v ; 

where is the Riemann zeta function of argument s. The mass-dependent coefficient 

-1 rl o,2n _ „l,2n _ „,2 , 

du I 

10 Jo 



A 2 (1/x) - / du I dv u2{l _ v2)+Ax2{l _ u) , W 



where x = mi/ m^ and mi is the mass of the virtual lepton in the vacuum polarization 
subgraph, was also computed in the late 1950s [19] for mi = m e and neglecting terms of 
0{m e /m^). Its exact expression was calculated in 1966 [20]. Actually, the full analytic 
result of [20] can be greatly simplified by taking advantage of the properties of the 
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dilogarithm Li 2 (z) = — J'q dt ln(l — t)/t. As a result of this simplification I obtain 



4 4, (l/s) = -— - — + x 2 (4 + 31nx)+x 4 



2 In x In ( x ) — Li2(x 2 ) 

3 \x 



+ |(l - 5x 2 ) y - lnx In (jq^) " Li 2 (x) + Li 2 (-x) 



(5) 



Note that this simple formula, contrary to the one in ref. [20], can be directly used 
both for < x < 1 (the case of the electron loop) and for x > 1 (tau loop). In the 
latter case, the imaginary parts developed by the dilogarithms Li 2 (x) and Li 2 (x 2 ) are 
exactly canceled by the corresponding ones arising from the logarithms. For x = 1 
(muon loop), eq. (5) gives 

4 4) (1) = H9/36 - ^ 2 /3; of course, this contribution is 
already part of in eq. (3). Evaluation of eq. (5) with the latest recommended 
values for the muon-electron mass ratio m M /m e = 206.768 2838(54) [21], and the ratio 
of = 105.658 3692(94) MeV [21] and m T = 1776.99(29) MeV [22] yields 

A ( 2 4 \m^/m e ) = 1.094 258 3111 (84) (6) 
A { 2 ] (m^/m T ) = 0.000078 064(25), (7) 

where the standard uncertainties are only caused by the measurement uncertainties 
of the lepton mass ratios. Eqs. (6) and (7) provide the first re-evaluation of these 
coefficients with the recently updated CODATA and PDG mass ratios of refs. [21,22]. 
These new values differ visibly from older ones (see refs. [8, 10]) based on previous 
measurements of the mass ratios, but the change induces only a negligible shift in 
the total QED prediction. Note that the r contribution in eq. (7) provides a ~ 42 x 
10~ n contribution to a^ ED . As there are no two- loop diagrams containing both virtual 

electrons and taus, i^fm^/me, m^/m T ) = 0. Adding up eqs. (3), (6) and (7) I get 
the new two-loop QED coefficient 

C 2 = + A { f{ mfx /m e ) + A^im^/mr) = 0.765 857410 (27). (8) 

The uncertainties in A^\mn/m e ) and A^' \m^/m T ) have been added in quadrature. 
The resulting error bC 2 = 2.7 x 10~ 8 leads to a tiny 0.01 x 10~ n uncertainty in a® ED . 



2.3 Three-loop Contribution 

More than one hundred diagrams are involved in the evaluation of the three-loop (sixth- 
order) QED contribution. Their analytic computation required approximately three 
decades, ending in the late 1990s. 

The coefficient A^ arises from 72 diagrams. Its calculation in closed analytic form 
is mainly due to Remiddi and his collaborators [23,24], who completed it in 1996 with 
the evaluation of the last class of diagrams, the non-planar "triple cross" topologies 
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Figure 2: The QED diagrams contributing to the muon g—2 in order a 2 . The mirror 
reflections (not shown) of the third and fourth diagrams must be included as well. 



(see, for example, fig. 3 A) [24]. The result reads: 



l 72 sv ; 24 sw 3 



° 4 + ~4 ( ln2 2 ~ 7I ' 2 ) ln2 2 



239 4 

-7T 



2160 



139 w . 298 „, 17101 9 28259 

+ C 3) vr 2 ln 2 + vr 2 + = 1.181 241 4566 . . . , (9 

18 sw 9 810 5184 ' v ' 

where a 4 = Y,n=i l/(2"n 4 ) = Li 4 (l/2) = 0.517479061 674 . . .. 

The calculation of the exact expression for the coefficient A 2 (m/M) for arbitrary 
values of the mass ratio (m/M) was completed in 1993 by Laporta and Remiddi [25,26] 
(earlier works include refs. [27,28]). For our analysis, m = m^, and M = m e or m T . 

This coefficient can be further split into two parts: the first one, A { 2 6) (m/M, vp), receives 
contributions from 36 diagrams containing electron or r vacuum polarization loops (see, 
for example, fig. 3 B) [25], whereas the second one, A 2 (m/M, lbl), is due to 12 light-by- 
light scattering diagrams with electron or r loops (like the graph of fig. 3 C) [26] . The 
exact expression for A 2 (m/M) in closed analytic form is quite complicated, containing 
hundreds of polylogarithmic functions up to fifth degree (for the light-by-light diagrams) 
and complex arguments (for the vacuum polarization contribution). It also includes 
harmonic polylogarithms [29]. As it is very lengthy, it was not listed in the original 
papers [25, 26] which provided, however, useful series expansions in the mass ratio 
(m/M) for the cases of physical relevance. Using the exact expressions in closed analytic 
form kindly provided to me by the authors, and the latest values for the mass ratios 
mentioned above, I obtain the following values 

^ 6) (m M /m e ,vp) = 1.920455130(33), (10) 
^ 6) (m M /m e ,lbl) = 20.94792489(16), (11) 
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A^(m M /m T ,vp) = -0.001782 33(48), (12) 

^> 6) (m M /m r ,lbl) = 0.002142 83(69). (13) 

The sums of eqs. (10)-(11) and eqs. (12)-(13) are 

Af\m^/m e ) = 22.868 380 02(20), (14) 

A^im^/nir) = 0.000360 51(21); (15) 



to determine the uncertainties, the correlation of the addends has been taken into 
account. Eqs. (10)-(15) provide the first re-evaluation of these coefficients with the 
recently updated CODATA and PDG mass ratios of refs. [21,22]. These new values 
differ visibly from older ones (see refs. [8, 10]) based on previous measurements of the 
mass ratios, but the change induces only a negligible shift in the total QED prediction. 
Note the large contribution from the electron light-by-light diagrams, eq. (11) - its 
leading term is (2/3)7r 2 ln(m M /m e ) [30]. More generally, it was shown in [31] that 
the 0(at 2n+1 ) contribution to a® ED , from diagrams in which the electron light-by- light 
subgraph is connected with 2n+l photons to the muon, contains a large TT 2n ln(m^/m e ) 
term with a coefficient of 0(1). 

The analytic calculation of the three-loop diagrams with both electron and r loop 
insertions in the photon propagator (see fig. 3 D) became available in 1999 [32]. This 
analytic result yields the numerical value 

4 6) (m M /m e ,m M /m T ) = 0.000 52766(17), (16) 

providing a small 0.7 x 10~ n contribution to a® ED . The error, 1.7 x 10~ 7 , is caused by 
the uncertainty of the ratio m M /m r . Combining the three- loop results presented above, 
I obtain the new sixth-order QED coefficient 

C 3 = 4 6) +4 6) (-M/-e)+4 6) K/^) + 4 6) (-M/-e,m /1 /rn r) 

= 24.050 509 64(43). (17) 

The error 5C 3 = 4.3 x 10~ 7 , due to the measurement uncertainties of the lepton masses, 
has been determined considering the correlation of the addends. It induces a negligible 
O(10~ 14 ) uncertainty in a® ED . 

In parallel to these analytic results, numerical methods were also developed, mainly 
by Kinoshita and his collaborators, for the evaluation of the full set of three-loop 
diagrams [7,33]. 

2.4 Four-loop Contribution 

More than one thousand diagrams enter the evaluation of the four-loop QED contribu- 
tion to a^. As only few of them are known analytically [34], this eighth-order term has 
thus far been evaluated only numerically. This formidable task was first accomplished 
by Kinoshita and his collaborators in the early 1980s [35,36]. Since then, they made a 
major effort to continuously improve this result [7,9,16,37-39], also benefiting from fast 
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Figure 3: Examples of QED diagrams contributing to the muon g— 2 in order a 3 . A, a 
"triple-cross" diagram. B (C), sixth-order muon vertex obtained by insertion of an electron 
or r vacuum polarization (light-by-light) subdiagram. D, graph with e and r loops in the 
photon propagator. 



advances in computing power. The latest analysis appeared in ref. [40]. One should 
realize that this eighth-order QED contribution, being about six times larger than the 
present experimental uncertainty of an, is crucial for the comparison between the SM 
prediction of and its experimental determination. 

There are 891 four-loop diagrams contributing to the mass-independent coefficient 

(S) (8) 

A\ . Its latest published value is A\ = —1.7502 (384) [39], where the error is caused 
by the numerical procedure. This coefficient has undergone a small revision in ref. [39]. 
In September 2004 Kinoshita reported a new preliminary updated result [41], 

Af ] = -1.7093(42). (18) 



Note the small shift in the central value and the significant reduction of the numerical 

(8) 

uncertainty of this new result. I will adopt it for the value of A\ . The latest value of 

(8) 

the coefficient (m^/rrie), arising from 469 diagrams, is [40] 



A^\m^/m e ) = 132.6823(72). (19) 

This value is significantly higher than the older one, 127.50 (41) [9] (its precision is 
impressively higher too) shifting up the value of a® ED by a non-negligible ~15 x 10~ n . 
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This difference is partly accounted for by the correction of a program error described 
in ref. [39], but is mostly due to the fact that the computation of the older value 
suffered from insufficient numerical precision. The term A 2 {m^/m T ) has been roughly 
estimated to give an O(10" 13 ) contribution to a® ED - it can be safely ignored for 
now [40]. The numerical evaluation of the 102 diagrams containing both electron and 
r loop insertions yields the three-mass coefficient [40] 

^ 8) (m M /m e ,m M /m r ) = 0.037594(83), (20) 

which provides a small O(10 -12 ) contribution to a® ED . Adding up the four-loop results 
described above, we obtain the eighth-order QED coefficient 

C 4 ~ 4 8) +A ( i\m lx /m e ) +Af\m lx /m e ,m lx /m T ) = 131.011(8). (21) 

(8) 

Note that this expression does not contain the term A 2 (m^/m T ), which has been 
roughly estimated to be of the same order of magnitude of the uncertainty on the r.h.s. 
of eq. (21). However, this uncertainty, 0.008, causes only a tiny 0.02 x 10~ n error in 



2.5 Five-loop Contribution 

The evaluation of the five-loop QED contribution is in progress [41]. The existing esti- 
mates are mainly based on the experience accumulated computing the sixth- and eighth- 
order terms, and include only specific contributions enhanced by powers of ln(m M /m e ) 
times powers of tt. The first estimate, C5 = 570(140), provided by Kinoshita and 
collaborators in 1990 [16], considered the contribution of graphs containing an electron 
light-by-light subdiagram with one-loop vacuum polarization insertions. A few other 
predictions for C5 exist, and classes of diagrams were computed or estimated with vari- 
ous methods [31,42-48]. In September 2004 Kinoshita reported a new very preliminary 
result [41], 

C 5 ~ 4 10) K/™e) = 677(40), (22) 

(9080 diagrams contribute to A^ \m^/m e )l) corresponding to a 4.6(0.3) x 10~ n 
contribution to a® ED . This is the value of C5 I will employ. The uncertainty in this 
new estimate of the tenth-order term (0.3 x 10~ n ) no longer dominates the error of 
the total QED prediction (see next section). Efforts to improve upon the evaluation of 
C5 are presently being pursued by Kinoshita and Nio. 



2.6 The Numerical Value of a^ ED 

Adding up all the above contributions and using the latest CODATA recommended 
value for the fine-structure constant [21], known to 3.3 ppb, 

a' 1 = 137.035 99911(46), (23) 
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I obtain the following value for the QED contribution to the muon g— 2: 

a QED = 116 584 718.8 (0.3) (0.4) x 10" 11 . (24) 

The first error is due to the uncertainties of the 0(a 2 ), 0(a 4 ) and 0(a 5 ) terms, and is 
strongly dominated by the last of them. (The uncertainty of the 0(a 3 ) term is negli- 
gible.) The second error is caused by the 3.3 ppb uncertainty of the fine-structure con- 
stant a. When combined in quadrature, these uncertainties yield 5a® ED = 0.5 x 10~ n . 
The value of a® ED in eq. (24) is close to that presented by Kinoshita in [41], a® ED = 
116 584 717.9 (0.3) (0.9) x 10~ n , and has a smaller error. This latter result was in fact 
derived using the value of a determined from atom interferometry measurements [49], 
a -1 = 137.036 000 3(10) (7.3 ppb), which has a larger uncertainty than the latest 
CODATA value employed for eq. (24). 



3 The Electroweak Contribution 

The electroweak (EW) contribution to the anomalous magnetic moment of the muon is 
suppressed by a factor (m M /M w ) 2 with respect to the QED effects. The one- loop part 
was computed in 1972 by several authors [50]. Back then, the experimental uncertainty 
of was one or two orders of magnitude larger than this one-loop contribution. Today 
it's less than one-third as large. 



3.1 One-loop Contribution 



The analytic expression for the one-loop EW contribution to a M , due to the diagrams 
in fig. 4, reads 



a EW (l loop) 



24v / 2~vr 2 



1 + i (l - 4sm 2 ^ w ) 2 + O 



m. 



(25) 



where = 1.16637(1) x 10~ 5 GeV -2 is the Fermi coupling constant, M z , M w and 
M H are the masses of the Z, W and Higgs bosons, and 9 W is the weak mixing angle. 
Closed analytic expressions for a EW (l loop) taking exactly into account the m 2 /M| 
dependence (B = Z, W, Higgs, or other hypothetical bosons) can be found in refs. [51]. 
Following [52], I employ for sin 2 ^ the on-shell definition sin 2 ^ = 1 - M^/Mj [53], 
where M z = 91.1875(21) GeV and M w is the theoretical SM prediction of the W mass. 
The latter can be easily derived from the simple analytic formulae of ref. [54], 



M w = 



80.4077 - 0.05738 In 



(on-shell scheme II with Aa[ 5) = 0.02761 (36), a s (M z ) = 0.118 (2) and M top = 178.0 (4.3) 
GeV [55]), leading to M w = 80.383 GeV for M H = 150 GeV, compared with the direct 
experimental value M w = 80.425 (38) GeV [22], which corresponds to a small M H [56]. 
For M H = 150 GeV, eq. (25) thus gives 



100 GeV 



0.00892 In 2 



100 GeV 



GeV, 



(26) 



a 



r {l loop) = 194.8 x 10" 



(27) 
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but this value encompasses the predictions derived from a wide range of values of M H 
varying from 114.4 GeV, the current lower bound at 95% confidence level [57], up to a 
few hundred GeV. 




Figure 4: One-loop electroweak contributions to a^. The diagram with a W and a Goldstone 
boson ((f>) must be counted twice. The diagrams with the Higgs boson loop and with two 
Goldstone boson couplings to the muon are suppressed by a factor m^/Mj w H and are not 
drawn. 



The contribution of the Higgs diagram alone, part of the 0(m 2 /M 2 WH ) terms of 



eq. (25), is [7,51] 



a.. • (1 loop) = — 7=r-^ 



logR H 7 +Q 



(28) 



where R H = M 2 /m 2 . Given the current lower bound M H > 114.4 GeV (95% CL), 
a^ w ' H (l loop) is smaller than 3 x 10~ 14 and can be safely neglected. 

3.2 Higher-order Contributions 

The two-loop EW contribution to was computed in 1995 by Czarnecki, Krause and 
Marciano [58,59]. This remarkable calculation, probably the first (and still one of the 
very few) complete two-loop electroweak computation, leads to a significant reduction 
of the one-loop prediction. Naively one would expect the two-loop EW contribution 
a^ w (2 loop) to be of order (a/ir) x a^ w (l loop), and thus negligible, but this turns out 
not to be so. As first noticed in 1992 [60], a^ w (2 loop) is actually quite substantial 
because of the appearance of terms enhanced by a factor of ln(M z ^ w /rrif), where mj 
is a fermion mass scale much smaller than M w . 

The two-loop contributions to a® w can be divided into fermionic and bosonic 
parts; the former includes all two-loop EW corrections containing closed fermion loops, 
whereas all other contributions are grouped into the latter. The full two-loop calcu- 
lation involves 1678 diagrams in the linear 't Hooft-Feynman gauge [61]. As a check, 
the authors of [58, 59] employed both this gauge and a nonlinear one in which the 
vertex of the photon, the W and the unphysical charged scalar vanishes. Their result 
for M H = 150 GeV (obtained in the approximation M H 3> M WiZ computing the first 
two terms in the expansion in M^ z /M 2 ) was a^ w (2 loop) = -42.3(2.0)(1.8) x 10~ u , 
leading to a significant reduction of a„ w . The first error is meant to roughly reflect low 
momentum hadronic uncertainties (more below), whereas the second allows for a range 
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of M H values from 114 GeV to about 250 GeV. Note that the contribution from 7— Z 
mixing diagrams is not included in this result: as it is suppressed by (1 — 4sin 2 # H/ ) ~ 0.1 
for quarks and (1 — 4sin 2 # xv ) 2 for leptons, it was neglected in this early calculation. 
It was later studied in ref. [52]: together with small contributions proportional to 
(1 — 4sin 2 # w ,)(m 2 /M{ 2 7 ) induced by the renormalization of sm 2 9 w , it shifts the above 
value of a^ w down by a tiny 0.4 x 10 -11 . 

The hadronic uncertainties, above estimated to be ~ 2 x 10~ n , arise from two types 
of two-loop diagrams: hadronic photon-Z mixing, and quark triangle loops with the 
external photon, a virtual photon and a Z attached to them (see fig. 5). The tiny 

Hadrons 




Figure 5: Hadronic loops in two-loop EW contributions. 



hadronic 7— Z mixing terms can be evaluated either in the free quark approximation or 
via a dispersion relation using data from e + e~ annihilation into hadrons; the difference 
was shown to be numerically insignificant [52] . The contribution from the second type of 
diagrams (the quark triangle ones), calculated in [58] in the free quark approximation, 
is numerically more important. The question of how to treat properly the contribution 
of the light quarks was originally addressed in ref. [62] within a low-energy effective 
field theory approach and was further investigated in the detailed analyses of refs. [52, 
63,64]. These refinements significantly improved the reliability of the fermionic part 
of a™ (two loop) and increased it by 2 x 10~ n relative to the free quark calculation, 
leading, for M H = 150 GeV, to [52] 

al w = 154(1)(2) x 10" 11 , (29) 

where the first error corresponds to hadronic loop uncertainties and the second to 
an allowed Higgs mass range of 114 GeV < M H < 250 GeV, the current top mass 
uncertainty 1 and unknown three-loop effects. 

The leading-logarithm three-loop contribution to a^ w was first studied via a renor- 
malization group analysis in ref. [65]. Such an analysis was revisited and refined in [52], 
where this contribution was found to be extremely small (indeed, consistent with zero 
to a level of accuracy of 10~ 12 ). An uncertainty of 0.2 x 10 -11 , included in eq. (29), has 

1 Indeed, although the result in eq. (29) was computed for ra t = 174.3 GeV, I checked that the shift 
induced in a^ n by the latest experimental value m t = 178.0 (4.3) GeV [55] is well within the quoted 
error. 
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been conservatively assigned to a^ w for uncalculated three-loop nonleading-logarithm 
terms. 

Lastly, I would like to point out that until recently only one evaluation existed 
of the two-loop bosonic part of a® w , ie, ref. [59]. The recent calculation of ref. [66], 
performed without the approximation of large Higgs mass previously employed, agrees 
with the result of [59] . Work is also in progress for an independent recalculation based 
on the numerical methods of refs. [67]. 

4 The Hadronic Contribution 

In this section I will analyze the contribution to the muon g— 2 arising from QED dia- 
grams involving hadrons. Hadronic effects in (two-loop) EW contributions are already 
included in a^ w (see the previous section). 

4.1 Leading-order Hadronic Contribution 

The leading hadronic contribution to the muon g— 2, a" LO , is due to the hadronic vac- 
uum polarization correction to the internal photon propagator of the one-loop diagram 
(see diagram in fig. 6). The evaluation of this 0(a 2 ) diagram involves long-distance 



Hadrons 




Figure 6: Leading hadronic contribution to a^. 

QCD for which perturbation theory cannot be employed. However, using analyticity 
and unitarity (the optical theorem) , Bouchiat and Michel [68] showed long ago that this 
contribution can be computed from hadronic e + e~ annihilation data via the dispersion 
integral [68,69] 

where a^(s) is the experimental total cross section for e + e~ annihilation into any 
hadronic state, with extraneous QED radiative corrections subtracted off (more later), 
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and R(s) is the ratio of cr(°\s) and the high-energy limit of the Born cross section 
for /i-pair production: R(s) = (s)/(4ira 2 /3s). The kernel K(s) is the well-known 
function 



Jo 



X 2 [L — X 



K{s) = / dx _ ■ (31) 



+ (1 — x)s/m 



(see ref. [70] for some of its explicit representations and their suitability for numerical 
evaluations). It decreases monotonically for increasing s, and for large s it behaves as 
m 2 /3s to a good approximation. For this reason the low-energy region of the dispersive 
integral is enhanced by ~ 1/s 2 . About 91% of the total contribution to a" LO is accumu- 
lated at center-of-mass energies y/s below 1.8 GeV and 73% of a" LO is covered by the 
two-pion final state which is dominated by the p(770) resonance [71]. Exclusive low- 
energy e + e~ cross sections have been mainly measured by experiments running at e + e~ 
colliders in Novosibirsk (OLYA, TOF, ND, CMD, CMD-2, SND) and Orsay (M3N, DM1, 
DM2), while at higher energies the total cross section ratio R(s) has been measured 
inclusively by the experiments 772, MARK I, DELCO, DASP, PLUTO, LENA, Crystal 
Ball, MD-1, CELLO, JADE, MARK-J, TASSO, CLEO, CUSB, MAC, and BES. Perturbative 
QCD becomes applicable at higher loop momenta, so that at some energy scale one can 
switch from data to QCD [72-74]. 

Detailed evaluations of the dispersive integral in eq. (30) have been carried out 
by several authors [70,71,73-93]. A prominent role among all data sets is played by 
the precise measurements by the CMD-2 detector at the VEPP-2M collider in Novosi- 
birsk [94-96] of the cross section for e + e~ — ► tt + tt" at values of yfs between 0.61 and 
0.96 GeV (ie, s G [0.37, 0.93] GeV 2 ). The quoted systematic error of these data is 
0.6% [96], dominated by the uncertainties in the radiative corrections (0.4%). In July 
2004, also the KLOE experiment at the DA<I>NE collider in Frascati presented the final 
analysis [97] of the 2001 data for the precise measurement of a(e + e~ — ► ir + ir~) via the 
radiative return method [98] from the (f> resonance. In this case the machine is operating 
at a fixed center-of-mass energy W ~ 1.02 GeV, the mass of the <f> meson, and initial- 
state radiation is used to reduce the invariant mass of the 7r + 7r~ system. In [97] the 
cross section a(e + e~ — ► ir + ir~) was extracted for the range s € [0.35,0.95] GeV 2 with 
a systematic error of 1.3% (0.9% experimental and 0.9% theoretical) and a negligible 
statistical one. The study of the e + e~ — > it + it~ process via the initial-state radiation 
method is also in progress at the BABAR detector at the PEP-II collider in SLAC [99]. 
This analysis will be important to further assess the consistency of the e + e~ data. The 
BABAR collaboration has already presented data for the it + it~it final state [100], and 
preliminary ones for the process e + e~ — > 2ir + 2ir~ [99]. On the theoretical side, the 
properties of analyticity, unitarity and chiral symmetry provide strong constraints for 
the pion form factor F n (s) in the low-energy region [92, 101-105]. They can lead to 
further improvements. Perhaps, also lattice QCD computations of a" LO , although not 
yet competitive with the precise results of the dispersive method, may eventually rival 
that precision [106]. 

The hadronic contribution a" LO is of order 7000 x 10 -11 . Of course, this is a 
small fraction of the total SM prediction for a M , but is very large compared with the 
current experimental uncertainty Sa^ xp = 60 x 10 -11 . Indeed, as 5a^ xp is less than 
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one percent of a^ LO , precision analyses of this hadronic term as well as full treatment 
of its higher-order corrections are clearly warranted. Normally, the "bare" cross section 
<t(°) (s) is used in the evaluation of the dispersive integral and the higher-order hadronic 
corrections (see sec. 4.2) are addressed separately. But what does "bare" really mean? 

The extraction of a^(s) from the observed hadronic cross section cr(s) requires the 
subtraction of several radiative corrections (RC) which, at the level of precision we are 
aiming at, have a substantial impact on the result. To start with, RC must be applied 
to the luminosity determination, which is based on large-angle Bhabha scattering and 
muon-pair production in low-energy experiments, and small-angle Bhabha scattering 
at high energies. The first step to derive a^°\s) consists then in subtracting the initial- 
state radiative (ISR) corrections (virtual and real, described by pure QED) from a(s). 
The resulting cross section still contains the effects of the photon vacuum polarization 
corrections (VP), which can be simply undressed by multiplying it by a 2 /a(s) 2 , where 
a(s) is the effective running coupling (obviously depending on nonperturbative contri- 
butions itself). The problem with data from old experiments is that it's difficult to find 
out if (and which of) these corrections have been included (sec rcf. [86]). The latest 
analysis from CMD-2 [96] is explicitly corrected for both ISR and VP (leptonic as well 
as hadronic) effects, whereas the preliminary data [94] of the same experiment were 
only corrected for ISR. For a thorough analysis of these problems, I refer the reader 
to [86,91,107] and references therein. 

All hadronic final states should be incorporated in the hadronic contribution to 
the muon g— 2, in particular final states including photons. These final-state radiation 
(FSR) effects, although of higher order (a 3 ), are normally included in the leading-order 
hadronic contribution a^ L °. I will stick to this time-honored convention. The precise 
CMD-2 data for the cross section e + e~ — > 7r + 7r~ (quoted systematic error of 0.6% 
dominated by the uncertainties in the RC) are corrected for FSR effects using scalar 
QED. I find this worrisome. The following is done: their experimental analysis imposes 
cuts to isolate the two-pion final states. These cuts exclude a large fraction of the 
7r + 7r~7 states, in particular those where the photon is radiated off at a relatively large 
angle [12]. The fraction left is then removed using the Monte Carlo simulation based 
on point-like pions. Finally, the full FSR contribution is added back using an analytic 
expression computed in scalar QED for point-like pions [108], shifting up the value of 
a* LO by ~ 50x 10" 11 [12,84,92,107,109,110]. (This full scalar-QED FSR contribution is 
also added to older 7r + 7r~ data.) This procedure is less than perfect, as it introduces a 
model dependence which could be avoided by a direct measurement of the cross section 
into hadronic states inclusive of photons. Any calculation that invokes scalar QED 
probably falls short of what is needed. 

The 2001 final analysis [95] of the precise CMD-2 tt + tt~ data taken in 1994-95 sub- 
stantially differed from the preliminary one [94] released two years earlier (based on the 
same data sample). The difference mostly consisted in the treatment of RC, resulting in 
a reduction of the cross section by about 1% below the p peak and 5% above. A second 
significant change occurred during the summer of 2003, when the CMD-2 collaboration 
discovered an error in the Monte Carlo program for Bhabha scattering that was used 
to determine the luminosity [96]. As a result, the luminosity was overestimated by 
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2-3%, depending on energy. (Another problem was found in the RC for /i-pairs pro- 
duction.) Overall, the pion-pair cross section increased by 2.1-3.8% in the measured 
energy range [71], a non-negligible shift. The 2004 results of the KLOE collaboration, 
obtained via the radiative return method from the <f> resonance, are in fair agreement 
with the latest energy scan data from CMD-2 [93,96,97]. Here I will only report the 
evaluations of the dispersive integral in eq. (30) based on the latest CMD-2 reanalysis, 
as it supersedes all earlier ones. These evaluations are in very good agreement: 2 

[93] a« LO = 6934(53)^(35)^ x HT 11 , (32) 

[88] af[ LO = 6948(86) x HT 11 , (33) 

[89] al h ° = 6934(92) x HT 11 , (34) 

[91] a« LO = 6924(59) e:cp (24) rad x HT 11 , (35) 

[92] a« LO = 6944(48)^(10)^ x HT 11 . (36) 

The preliminary result in eq. (32) already includes KLOE's 2004 data analysis and up- 
dates the one of ref. [71], shifting it down by 29 x 10 -11 ; two thirds of this shift are due 
to the inclusion of KLOE's data. The preliminary new result in eq. (34) updates the 
value a^ LO = 6996 (85) e:r p(19) ra d(20)p roc x 10 -11 previously obtained by the same au- 
thors [90] . Their central value decreased because of an improvement of their integration 
procedure. 

The authors of ref. [81] pioneered the idea of using vector spectral functions derived 
from the study of hadronic r decays [111] to improve the evaluation of the dispersive 
integral in eq. (30). Indeed, assuming isospin invariance to hold, the isovector part 
of the cross section for e + e~ — > hadrons can be calculated via the Conserved Vector 
Current (CVC) relations from r-decay spectra. An updated analysis is presented in [71], 
where r spectral functions are obtained from the results of ALEPH [112], CLEO [113] 
and OPAL [114], and isospin-breaking corrections are applied [115-117]. In this r-based 
evaluation, the 2ir and the two 47r channels are taken from r data up to 1.6 GeV and 
complemented by e + e~ data above (the QCD prediction for R(s) is employed above 
5 GeV). Note that r decay experiments measure decay rates which are inclusive with 
respect to radiative photons. Their result is 

[71] a» LO = 7110 (50)^(8)^(28)5^(2) x 10 _1 \ (37) 

where the quoted uncertainties are experimental, missing radiative corrections to some 
e + e~ data, and isospin violation. This value must be compared with their e + e~-based 
determination in eq. (32). Also the analysis of [92] includes information from r decay. 
They obtain 

[92] a» LO = 7027(47) ea;p (10) ra( i x 10- 11 , (38) 

to be compared with their determination in eq. (36). 

2 I have translated the results of ref. [92] into the notation of the present article. 
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Although the latest CMD-2 e + e~ — ► ir + ir~ data are consistent with r data for the 
energy region below 850 MeV, there is an unexplained discrepancy for larger energies. 
This is clearly visible in fig. 7, from ref. [93], where the relative comparison of the 7r + 7r~ 
spectral functions from e + e~ and isospin-breaking-corrected r data is illustrated. The 
same figure also shows the tt + tt~ spectral functions derived from KLOE's 2004 e + e~ 
analysis. They are in fair agreement with those of CMD-2 and confirm the discrepancy 
with the r data. 
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Figure 7: Relative comparison of the it + it~ spectral functions from e + e~ and isospin- 
breaking-corrected r data, expressed as a ratio to the r spectral functions. The band shows 
the uncertainty of the latter. This figure is from ref. [93]. 



Among the possible causes of this discrepancy, one may wonder about inconsisten- 
cies in the e + e~ data, in the r data, or in the isospin-breaking corrections applied to 
the r spectral functions. Given the good consistency of the ALEPH and CLEO data 
sets, and the confirmation by KLOE of the trend exhibited by other e + e~ data, fur- 
ther careful investigations of the isospin-violating effects are clearly warranted - see 
the interesting studies in [92,93,103,118-120], in particular the discussion of the pos- 
sible difference between the masses and the widths of neutral and charged p-mesons. 
Until we reach a better understanding of this problem, it is probably safer to discard 
information from r decays for the evaluation of a^ LO [93]. 
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4.2 Higher-order Hadronic Contributions 

We will now briefly discuss the 0(a 3 ) hadronic contribution to the muon g— 2, a" HO , 
which can be divided into two parts: 

<" = «r°(vp) + <- (lbl). (39) 

The first term is the 0(a 3 ) contribution of diagrams containing hadronic vacuum po- 
larization insertions, including, among others, those depicted in figs. 8 A and B. The 
second one is the light-by-light contribution, shown in fig. 8 C. Note that the 0(a 3 ) 
diagram in fig. 8 D has already been included in the leading-order hadronic contribu- 
tion a^ LO although, unsatisfactorily, using scalar QED (see discussion in sec. 4.1). In 
recent years, a^ HO (vp) was evaluated by Krause [121] and slightly updated in [81]. Its 
latest value is [91] 

< ff °(vp) = -97.9 (0.9) e!rp (0.3) ro d x 1(T U . (40) 

This result was obtained using the same hadronic e + e~ annihilation data described in 
sec. 4.1. It changes by about —3 x 10~ n if hadronic r-decay data (again, see sec. 4.1) 
are used instead [15]. 




Figure 8: Some of the higher-order hadronic diagrams contributing to a M . 



The hadronic light-by-light contribution changed sign already three times in its 
troubled life. Contrary to a^ HO (vp), it cannot be expressed in terms of experimental 
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observables determined from data and its evaluation therefore relies on purely theoreti- 
cal considerations. The estimate of the authors of [13, 122, 123], who uncovered in 2001 
a sign error in earlier evaluations, is 

[13] < H °(lbl) = + 80(40) x lO" 11 . (41) 

Earlier determinations now agree with this result [124,125]. Further studies include [126, 
127]. At the end of 2003 a higher value was reported in [128], 

[128] <"°(lbl) = +136(25) x HT 11 . (42) 

It was obtained by including short-distance QCD constraints previously overlooked. 
Further independent calculations would provide an important check of this result for 
a" HO (lbl), a contribution whose uncertainty may become the ultimate limitation of the 
SM prediction of the muon g— 2. 

5 The Standard Model Prediction vs. Measurement 

We now have all the ingredients to derive the SM prediction for a M : 

< M = af° + 4 W + a T° + < H °M + < H °( lb1 )- ( 4 3) 
For convenience, I collect here the values of each term from eqs. (24, 29, 32-38, 40-42): 
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The values I obtain for a^ M are shown in the first column of table 1. The values em- 
ployed for a" LO are indicated by the reference in the last column. I used the latest value 
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available for the hadronic light-by-light contribution a" HO (Vo\) = 136 (25) x 1CT 11 [128]. 
Errors were added in quadrature. 

The latest measurement of the anomalous magnetic moment of negative muons by 
the experiment E821 at Brookhaven is [5] 

a Exp = H6 592 140 (80) (30) x 10" 11 , (44) 

where the first uncertainty is statistical and the second is systematic. This result is 
in good agreement with the average of the measurements of the anomalous magnetic 
moment of positive muons [2-4, 129], as predicted by the CPT theorem [130]. The 
present world average experimental value is [5] 

a° xp = 116 592 080(60) x 10~ n (0.5 ppm). (45) 

The comparison of the SM results with the present experimental average in eq. (45) 
gives the discrepancies {a E \ xp — a„ M ) listed in the second column of table 1. The 
number of standard deviations, shown in the third column, spans a wide range from 
0.7 to 2.8. Somewhat higher discrepancies, shown in parentheses in the third column, 
are obtained if the hadronic light-by-light contribution a" HO (lbl) = 80 (40) x 10~ n [13] 
is used instead of af[ HO (M) = 136(25) x 10~ n [128], with the number of standard 
deviations spanning the range [1.3 — 3.2] instead of [0.7 — 2.8]. Note that the entries of 
the first row in table 1 are based on the preliminary result for a" LO of ref. [93], which 
already includes the recent data from KLOE and updates the one of ref. [71], shifting it 
down by 29 x 10~ n . As two thirds of this shift are due to the inclusion of the KLOE data, 
it is possible that eventually also the a" LO results of refs. [88,89,91,92] will undergo 
some decrease as a consequence of this inclusion, thus increasing the corresponding 
af l xp — a^ M discrepancies. 



al M x 10 11 


K XP ~ < M ) x 
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HLO Reference 


116591845 (69) 
116591859 (90) 
116591845 (95) 
116591835 (69) 
116591855 (55) 


235 (91) 
221 (108) 
235 (113) 
245 (91) 
225 (81) 


2.6 (3.0) 
2.1 (2.5) 
2.1 (2.5) 

2.7 (3.1) 

2.8 (3.2) 
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[92] (e+e-) 


116592018 (63) 
116591938 (54) 


62 (87) 
142 (81) 


0.7 (1.3) 
1.8 (2.3) 


[71] (r) 

[92] (e+e-,r) 



Table 1: SM predictions for compared with the current measured world average value. 
See text for details. 

6 Conclusions 

In the previous sections I presented an update and a review of the contributions to 
the SM prediction for the muon g— 2. What should we conclude from the wide spec- 
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trum of results obtained in sec. 5? The discrepancies in table 1 between recent SM 
predictions and the current world average experimental value range from 0.7 to 3.2 
standard deviations, according to the values used for the leading-order and light-by- 
light hadronic contributions. In particular, the contribution of the hadronic vacuum 
polarization depends on which of the two data sets, e + e~ collisions or r decays, are 
employed. 

This puzzling discrepancy between the ir + ir~ spectral functions from e + e~ and 
isospin-breaking-corrected r data could be caused by inconsistencies in the e + e~ data, 
in the r data, or in the isospin-breaking corrections applied to the latter. Given the 
fair agreement between the CMD-2 and KLOE e + e~ data, and the good consistency of 
the ALEPH and CLEO r spectral functions, it is clear that further careful investigations 
of the isospin violations are highly warranted. Indeed, the question remains whether 
all possible isospin-breaking effects have been properly taken into account. Until we 
reach a better understanding of this problem, it is probably safer to discard information 
from hadronic r decays [93]. (Of course, discarding r data information still leaves us 
with the problem of their discrepancy, a troublesome issue on its own, independent of 
the calculation of the muon g—2.) If e + e~~ annihilation data are used to evaluate the 
leading hadronic contribution, the SM prediction of the muon g—2 deviates from the 
present experimental value by 2-3 standard deviations. 

The measurement of the muon g—2 by the E821 experiment at the Brookhaven 
Alternating Gradient Synchrotron, with an impressive relative precision of 0.5 ppm, is 
still limited by statistical errors rather than systematic ones. A new experiment, E969, 
has been approved (but not yet funded) at Brookhaven in September 2004 [131]. Its 
goal would be to reduce the present experimental uncertainty by a factor of 2.5 to about 
0.2 ppm (±23 x 10 -11 ). A letter of intent for an even more precise g—2 experiment was 
submitted to J-PARC with the proposal to reach a precision below 0.1 ppm [132]. While 
the theoretical predictions for the QED and EW contributions appear to be ready to 
rival these precisions, much effort will be needed in the hadronic sector to test a^ M at an 
accuracy comparable to the experimental one. Such an effort is certainly well motivated 
by the excellent opportunity the muon g—2 is providing us to unveil or constrain "new 
physics" effects. 
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